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Ph' Abstract 

^ . Recently measured VHE gamnia ray spectra from supernovae remnants (SNRs) are 

c/3 , best fit by power laws with an exponential cutoff. But this feature does not occur at 

^ ' the 'knee-equivalent' energy at which VHE gamma ray spectra are expected to reflect 
the 'knee' in the otherwise featureless cosmic ray spectrum. In this article, the VHE 

^ ! gamma ray feature is explained as a consequence of a recently deduced quantum effect. 

X/-^ ' The effect distinguishes 'trajectory energy' from 'particle state energy' and requires 

00 ■ that the particle state energy depends strongly on gravitational potential at very high 

^y-^ . energies. Based on this effect and the observed CR spectrum, a tight two-parameter fit 

J.." I is obtained to the combined VHE gamma ray spectra of SNRs RX J1713. 7-3946 and 

^ ■ RX J0852. 0-4622 and the Galactic Center ridge. 
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1 INTRODUCTION 2 

1 Introduction 

It is widely held that galactic cosmic rays are mostly atomic nuclei accelerated in supernovae 
remnants, that some CRs collide with local ambient matter producing gamma rays, and that 
CRs, being charged, are deflected as they travel to Earth. Being uncharged, the gamma rays 
that arrive at Earth retain directional information pointing back to their sources. 

Since the galactic CR primaries have been deflected, the CR spectrum observed at Earth, 
see Fig. [H represents an amalgam of SNR sources. The spectrum follows a near-constant 
power law from 10^° eV to the 'knee' at 3 x 10^^ eV, where the spectral index of the power 
law steepens and is again constant to the 'ankle' at 10^*^ eV. The 'ultrahigh energy' CRs do 
not concern us here. The most energetic VHE gamma rays thus far observed have energies 
10^^ eV or so, implying that the relevant CR energies are below about 10^^ eV. 

Considering only proton primaries, CR protons at the SNR are expected to produce VHE 
gamma rays when the protons collide with ambient matter, most likely protons. The p-p 
collisions create pions and etas that decay to produce the VHE gamma rays. The process is 
well understood and the only quantity that could provide the VHE gamma ray spectra with 
a distinct feature is the inelastic p-p cross section ainei- 

However the log of ainei is a slowly rising function of incident proton energy and is 
remarkably featureless, see Fig. [3^. Thus one expects the VHE gamma ray spectra from CR 
protons to be a power law up to some 'knee equivalent' gamma ray energy, and at higher 
energies to again be a power law but with a probably steeper index. In short, the VHE 
gamma ray spectrum should be similar to the CR spectrum and a typical VHE spectrum 
should feature a knee equivalent somewhere around 0.17 x 3 x 10^^ = 5.1 x 10^^ eV. 

The H.E.S.S. collaboration reports a cutoff in the spectrum for the supernova remnant 
RX J1713. 7-3946 (Aharonian et al. 2007a) at an energy of at most 1.8 x 10^^ eV, a factor 
of about 30 below the expected knee equivalent. Therefore, it is thought that the feature 
reported is not the knee equivalent (Plaga 2008a). One must consider other explanations. 
It could be that the VHE gamma rays are emitted by electrons in the 'lepton scenario' since 
electrons are also believed to be accelerated in supernovae remnants and electrons in certain 
conditions emit VHE gamma rays. For lepton scenario discussions, see Huang et al (2007), 
Katz & Waxman (2008), and Plaga (2008b). 

It is shown in this paper that the hadron scenario fits the VHE gamma ray feature at 
10^'^ eV if we apply the result of a recent calculation. The explanation distinguishes two 
different kinds of energy that are equal at low energies and diverge at very high energies. 
The 'particle state energy' depends on gravitational potential, with less energy at the SNR 
and more energy at Earth (Shurtleff 2008). A summary of the calculation is included in 
Sec. [2] for convenience. 

Then a massive particle such as a proton has two energies, one that can be obtained by 
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mass and a time-of-flight velocity, a 'trajectory energy'. The other is the energy dehvered 
in a colhsion due to the change in the particle's state. Call that the 'particle state energy'. 
The formula for the particle state energy E is 

E^ 
E = Eo-A<p^, (1) 

where < is the gravitational potential and Eg is the energy at null potential = 0. The 
gravitational field is assumed to be weak and the potential is the Newtonian potential 
made unitless by dividing by the square of the speed of light c. The mass of the particle is 
m. 

The Solar System moves in the gravitational field of the Galaxy. Very roughly, the Galaxy 
should contribute a potential (pcai due to Mcai = 10^^ solar masses in effect collected some 
R = 2.5 X 10^° m away from the Solar System at the Galactic Center. One has (pcai = 
—GMcai/Rc^ ^ — 6 X 10"'^. It follows that the effects of Eq. (^ would be significant at 
about Eq = mc^/2y/—(f)aai = 0.6 TeV, where mc^ ^ 0.001 TeV is the rest energy of a proton. 
Since 0.6 TeV is well within the range of proton energies available at the Tevatron, and no 
such effects are reported from the Tevatron, the highest energy proton accelerator as of this 
writing, it follows that the potential due to the Galaxy and any larger potentials must be 
eliminated from consideration if Eq. ([1]) is true for protons. 

To eliminate the effects of the Galaxy's gravitational field, one may choose to work in 
a comoving coordinate system in free fall with the Solar System, see for example Weinberg 
(1972). In General Relativity, a particle's trajectory can be described in any coordinate 
system. 

By making the particle state energy dependent on potential and requiring the potential 
to be calculated in a system in free fall with the Solar System, the description is valid in 
one coordinate system, the reference frame of the Solar System. This is allowed in General 
Relativity as long as the various quantities are properly transformed when the description 
is given in some other coordinate system. 

While the coordinate system is in free fall with the Solar System, thereby canceling the 
Galactic and larger scale potentials, the local gravitational potential due to the matter in 
the Solar System does not cancel. Thus, at the Earth's surface, the in ([1]) is the potential 
due to the Sun and the Earth, a value of about 00 ~ — 1 x 10~^, which becomes significant 
at an energy Eq ^ 5 TeV. This is beyond the proton energies available at proton accelerators 
at the time of this writing, but is within the reach of the LHC which has not begun high 
energy runs at the time of this writing. 

While the properties of electrons have no bearing on the explanation of the feature in 
the VHE gamma ray spectrum, one may note the energy at which the effects of ([T]) would 
become significant for electrons is Eq = meC^/2A/— 0^ = 2.5 GeV. Since the effects expected 
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by applying ([1]) have not been reported by accelerator experiments with electron energies 
many times the threshold 2.5 GeV, it follows that electrons do not obey ([T]), at least for 
potentials (p with magnitudes as large as the potential due to the Earth and Sun. This could 
mean the comoving coordinate system for electrons is in free fall with matter falling in the 
gravitational field due to the Earth and Sun (and due to the Galaxy and larger scales as 
well). Or it could mean that the assumptions that produce ([T]) do not apply to electrons. 
This could be another way that electrons differ from protons. 

One benefit of going to a system in free fall with the local matter is the resulting univer- 
sality possible for the various SNRs. If the potential at each SNR includes the contribution 
from the Galaxy, then Eq. ([1]) would have a different (ps for each supernova remnant. That 
would be a troubling problem. But if the potential (ps is the local potential in the system in 
free fall with the SNR, then the potentials at the various SNRs are likely to be nearly the 
same and Eq. ([1]) can meaningfully be applied to SNRs anywhere in the Galaxy. 

In this picture, to get a cosmic ray from the SNR to Earth, one must consider a number 
of overlapping patches of spacetime, each patch described by a coordinate system in free fall 
with the local matter at a scale comparable to objects like the Solar System or supernova 
remnants. 

With Eq. ([T]), a cosmic ray proton observed at Earth with energy Ep(^ in the potential 
00 at the Earth has a different energy Eps at the supernova in the potential (ps- In this 
notation, 'p' stands for proton and '©' for Earth and 'S" for supernova remnant. 

Relating the CR spectrum at Earth, dNp/dEp^, to the CR spectrum at a typical SNR, 
dNp/dEps, involves the derivative dEp^/dEps relating the energy at the SNR and the 
corresponding energy at Earth. It is suggested here that this energy dependence is part of 
the feature found by the H.E.S.S. collaboration in the VHE spectrum of SNR RX J1713.7- 
3946. 

Having obtained the expected CR spectrum at a typical SNR based on the CR spectrum 
at Earth, the expected VHE gamma ray spectrum at the SNR can be determined by the 
analysis by Kelner et al. (2006). Then the expected VHE gamma ray spectrum at Earth 
can be found if we know how it is related to the VHE gamma ray spectrum at the SNR. 

The problem is how to relate the VHE gamma ray spectrum at the SNR to the VHE 
gamma ray spectrum at Earth. Since the basic equation ([1]) applies to massive particles, one 
may entertain various arguments to get it to apply to massless gamma rays. 

First of all, note that at 10 TeV, say, a proton's rest energy is just 0.01% of its total 
energy. The proton moves very much like a massless photon of the same energy. So a 10 
TeV VHE gamma ray should behave much like a 10 TeV proton. Since the basis of the 
effect separates this 'trajectory energy' from the 'particle state energy', this argument may 
be spurious, even though it is intuitively appealing. 

Secondly, consider the process of neutral pion decay that produces many of the VHE 
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gamma rays in the hadronic scenario. A 10 TeV neutral pion decays to a 10 TeV VHE 
gamma ray plus a 4 keV X-ray. Within roundoff, the pion and the gamma ray each have an 
energy of 10 TeV. By time inversion, starting with a 10 TeV VHE gamma ray, an interaction 
with a 4 keV X-ray could produce a 10 TeV neutral pion. However unlikely the process, it 
could happen, so a 10 TeV VHE gamma ray is virtually a 10 TeV pion. Since the virtual 
transition from VHE gamma ray to neutral pion could occur anywhere from the SNR to 
Earth and a neutral pion traveling from the SNR to Earth would obey Eq. ([T]), the energy 
of the VHE gamma should obey Eq. ([1]) with m as the neutral pion mass 771^^0. 

Imagine placing a shuttered X-ray source at each kilometer along the path of a VHE 
gamma ray from an SNR to Earth. Furthermore let there be a limitless supply of identical 
VHE gamma rays ready to be sent from the SNR to Earth one by one. Let the n = VHE 
gamma ray travel from the SNR to Earth without removing any of the shutters. This VHE 
gamma ray reaches Earth unchanged by the presence of the X-ray sources. For the next VHE 
gamma ray, i.e. n = 1, release X-rays at the first kilometer at the SNR so that a neutral pion 
is created by inverse pion decay. The neutral pion has almost the same energy as the VHE 
gamma ray. For the nth VHE gamma ray release X-rays from just the one source at the nth 
kilometer making an nth neutral pion which has almost the same energy as the nth VHE 
gamma ray at the nth kilometer. Since the n = VHE gamma ray has the same energy as 
the nth VHE gamma ray just before the nth kilometer, and the nth VHE gamma ray has 
nearly the same energy as a neutral pion at the nth kilometer, and a neutral pion obeys ([1]), 
it follows that the VHE gamma ray energy should also obey Eq. ([1]) with m equal to the 
mass of a neutral pion. 

It should be clear that the result does not depend on there being any real X-rays along 
the path of the VHE gamma ray from the SNR to Earth. It should also be clear that none 
of this implies that a VHE gamma ray has mass, only that its energy could be transferred 
essentially undiminished to a neutral pion in some possibly imaginary circumstances. 

But there are other ways a 10 TeV VHE gamma ray could interact with ambient matter 
of negligible energy, perhaps making the VHE gamma ray virtually an electron or a positron 
or an eta meson, etc. Again, there need not be any real ambient matter, imaginary ambient 
matter suffices. Not knowing how to calculate the likelihood of independent imaginary 
interactions means not knowing what mass to associate with the VHE gamma ray. Thus we 
let that mass be a free parameter. 

By these arguments, we apply Eq. ([T]) to VHE gamma rays by letting the mass m = 
rriray staud for the average mass of the masses of massive particles that could be given the 
same energy as the VHE gamma ray by some possibly imaginary process. The mass rriray is 
a free parameter found by fitting the experimental VHE gamma ray spectra. 

Allowing ([1]) to apply to the VHE gamma rays allows us to use the expected VHE gamma 
ray fiux at the SNR, dN^/dE'-fs, to determine the expected VHE gamma ray fiux at Earth, 
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dN^/dEji^. The relationship involves the derivative dEjs/dE'jt^ which is energy dependent. 
As discussed previously, one also has the derivative dEp^/dEps from the CR spectrum. 
The product of the two factors, dEp^/dEps times dE'js/dE'y^, is an energy dependent 
quantity, see Fig. [3)d, that could describe the feature in the VHE gamma ray spectrum of 
the SNR RX J1713.7-3946 reported by the H.E.S.S. collaboration. 

2 Potential-Dependent Particle State Energy 

This section is included to indicate the origin of the basic formula. For details see Shurtleff 
(2008). 

The quantum field for a massive particle species can be expanded over the annihilation 
and creation operators that remove or add particle states to multiparticle states, see, for 
example, Weinberg (1995). Consider an annihilation field ?/'"'", 

'^^i^) =^ j d^V uUx, 'f)a„{;f) , (2) 

where the annihilation operator a^i^^ removes a state of the particle with (four-)momentum 



{"p^,p*} with total energy p* = +\Jm?' — jf . The ui„{x^ p^) are the 'coefficient functions.' 

To comply with special relativity, the quantum field ip'j^{x) transforms linearly under a 
Poincare transformation {A, 6}, where A is the result of successive rotations and boosts, a 
'Lorentz transformation', followed by a translation along the displacement b. Such transfor- 
mations preserve the scalar products of coordinate differences between events. Likewise the 
operator a^ilf) transforms linearly with a Poincare transformation. 

An essential step in obtaining the basic formula needed in this article depends on a char- 
acteristic of Poincare transformations, a property of translations that is not shared with 
rotations and boosts. Unlike with rotations and boosts, under a translation the coordinate 
differences themselves are invariant. Thus it is impossible to tell by keeping track of co- 
ordinate differences just what the displacement was. One displacement e could have been 
applied to the annihilation operator and a displacement b could have been applied to the 
quantum field. The displacements e and b may be independent of one another without any 
effect on the underlying coordinate differences. 

The annihilation and creation operators and the particle states must transform via uni- 
tary reps to satisfy fundamental quantum principles. The transformation {A, e} applied to 
the annihilation operators can be written as follows. 



f/(A, e)a^it)U-\A, e) = e^^^^^ ^ DiT'a.irp) , (3) 
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where D^^' is the transformation matrix and j is the spin of the particle. What is important 
here is the phase factor, e*^^*^, which is a unitary representation of {A, e}. 

With the transformation {A, e} apphed to the annihilation operators, the quantum field 
is required to transform as follows, 

UiA,b)^+ix)U-\A,b)=Y.D;^\A,b)^fiAx + b), (4) 

I 

where Ax + 6 results from transforming the coordinates and -D(A, b) is the nonunitary trans- 
formation matrix corresponding to the possibly reducible spin, say {A, B) © (C, D) © . . . , of 
the field %p^ . For a spin j = 1/2 Dirac particle the field transforms with spin (1/2, 0)©(0, 1/2). 

Thus one object, the quantum field iIj~^, is a linear combination of other objects, the 
operators ao-("p^), that transform differently from the first type of object. In such cases, the 
coefficients, here the uia{x, Jf), are constrained. Clebsch- Gordon coefficients can be derived 
from such constraints involving rotations in Euclidean 3-space (Shankar 1994). Here the 
transformations differ in part because one is unitary and the other nonunitary. Also the 
displacement for the operators e may differ from the displacement b for the quantum field. 
Deriving free quantum fields from special relativity in this way is presented in Chapter 5 of 
Weinberg (1995). 

In order to obtain the basic formula, consider the arbitrary displacement e. Since e is 
arbitrary, it may be useful to restrict it somewhat and make it depend on various relevant 
quantities. Let e depend on the field's Poincare transformation {A, b} and also depend on 
the event coordinates x. The constraints that determine the coefficient functions uia{x, Jf') 
also constrain the function e{A,x,b). 

Leaving the details to the references, one finds that the coefficient functions are deter- 
mined by (0), dSD, Q to be 

uUx, -f) = J- e^"-*''^ E DiiiL, x)uUO, ) , (5) 

\P I 

where the coefficient functions on the right are evaluated at the origin x = 0. The Lorentz 
transformation L transforms the rest momentum { ,m} to the momentum p^ = {'jf^,p^}- 
The quantity M = M{x) is a second rank tensor field that is arbitrary until the motion of 
the particle is considered. 

The displacement function e{A,x,b) is also constrained by Eq. ([2]), ([3]) and (jl]). One has 

e'^(A, X, b) = e^{A) - A'^[M{x)YX + [M{Ax + bX{Ax + bf , (6) 

where M{x) is the arbitrary second order tensor field that also appears in ([5]). We drop the 
first term, e^(A) = 0. The displacement e is equal to b when M = 1, i.e. when M is the 
delta function, the displacements of the field and operators are equal. 
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The classical motion of a particle follows paths of extreme quantum phase. The relevant 
phase is the phase of the quantum field, 6Q = p ■ M6x. One finds that 6Q is extreme when 
M6x is parallel to the momentum p. One may write this as p" = mM^dx^/dr where the 
quantity r turns out to be the proper time. Since the momentum is related to mass by 

?7a/3P"p^ = -m^ , (7) 

one finds immediately from p" = mM'^dx^ / dr that 

dx^^ dx" , , 

where the 'curved spacetime metric' g^^y is defined by 

^M^^^-z^M^Mf. (9) 

Further calculations show that the extreme path is the path from general relativity for the 
curved spacetime metric g^i,y (Shurtleff 2008). Since the tensor M is arbitrary, it follows 
that the metric g is arbitrary and the trajectory can be described in any coordinate system 
limited only by conditions of continuity and differentiability that are not considered here. 

By ([9]), the second rank tensor field M[x) is a kind of square root field of the curved 
spacetime metric g^y[x). In a weak gravitational field, to sufficient accuracy, in terms of the 
Newtonian potential 0, the metric gfj.u{x) is, see for example Adler et al. (1965), 

gf,u = diag{gxx, gxx, gxx, gtt) = diag(l - 20, 1 - 20, 1 - 20, -1 - 20) . (10) 

Then, by ([H]), one choice for M is 

M; = diag(M,, M,, M,, Mt) = diag(l - 0, 1 - 0, 1 - 0, 1 - 0, 1 + 0) , (11) 

where terms of second order in are dropped. In this way the classical motion determines 
the value of the field M{x). 

For the Newtonian potential at Earth we have 

= 0e + 0^ = -^r^-Y^^ + 0^, (12) 

where G is the universal gravitational constant, M^ and Mq are the masses of the Earth 
and Sun. Also, r^ is the radius of the Earth, c is the speed of light and 0oo is the potential 
far from the Solar System and far from any other massive object. Note that the 0s have 
been made unitless by dividing by c^. 
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In order to equate the potential (j) with cj)^ in the expression ( ITOl) for the curved spacetime 
metric g^j,^, it is necessary to make the potential 0oo vanish. To do this one can insist that the 
trajectory be described in a comoving coordinate system in free fall with the Solar System. 
In any such coordinate system, the time-time component of the curved spacetime metric 
is determined, gu = —1, (Weinberg 1972). In weak fields this component depends on the 
Newtonian potential gu = —1 — 2(j)oo- Thus in a comoving coordinate system in free fall with 
the Solar System and at a large distance from the Solar System, (poo = 0. Then the potential 
0OO of the Galactic gravitational field and more distant masses is canceled out, leaving the 
local gravitational potential to dominate near the Solar System. Thus Eq. (1 10 1) applies with 
= 0® in the comoving coordinate system in free fall with the Solar System. 

These considerations depend on the phase of the quantum field which, for a given mo- 
mentum p, is the phase p ■ Mx of the associated coefficient function uia-{x, Jf). What about 
the phase of the particle state with momentum p7 

First of all, note that the operator ao-(Ap) removes the same state no matter where the 
field ip~^{x) is defined, so the operator is independent of x. Since the operator does not 
communicate location information between the field and the particle states, one cannot be 
sure that the annihilation operator acr{'p') removes a state of the particle with momentum 
{"p*',p*} in the same reference frame as the quantum field is defined in. Therefore, one can 
allow the particle states to be defined in a reference frame with coordinates y while the field 
is defined in another frame at events with coordinates x. It follows that x = Xy, for some 
Lorentz transformation A. 

The transformation {A, e} changes the initial phase po-yo of the particle state with initial 
momentum p^ as follows: 

g«Po-J/o _^ ^ip-{y+f-) — ^ip-y^ip-(,-M+A^''-MA)\y^ip-MoSx j-j^gN 

where the last expression follows from (jS]). Also, p = ApQ, y = Ayo, Mq is the tensor in 
the initial reference frame and M is in the transformed frame with Mq = A~^MA and M = 
AMoA-\ 

For a particle moving in a weak gravitational field, one can treat M as constant over 
fairly large regions, large on the scale of the relevant portion of the quantum field. And on 
such a scale, the change of the momentum p is often small due to the weak gravitational 
force. Thus, in many cases, we can treat both M and p as constants on a scale much larger 
than the scale of the quantum field. 

Then the coefficient of y in the phase of the exponential in (fT3|) does not depend on y 
and we have a plane wave in a region where M and p can be considered constant. Defined 
as a constant times the rate of change of phase with distance, the momentum of the particle 
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state's plane wave is not p, but an 'effective momentum' p given in 

p.y = p.y + p. (-M + A-^MA)Xy = p{l - MX + A^^MAA) ■ y , 

from wliicli it follows that 

p = p{l-M\ + A-^MAX) (14) 

The two momenta p and p are equal when A = 1 or M = 1. 

Since the momentum of the particle state plane wave depends on the Lorentz transfor- 
mations A and A, it is convenient to work in a fixed reference frame. Appropriate transfor- 
mations can be applied as needed to go to other frames. A suitable fixed frame is provided 
by the distribution of Cosmic Microwave Background (CMB) radiation. The observed dipole 
anisotropy of the CMB implies the Solar System is moving at a speed of 370 km/s = 0.00123c 
(Fixsen et al. 1996, Lineweaver et al. 1996) with respect to the distribution of CMB radi- 
ation. We take the given frame to be the CMB reference frame and, by ignoring the 370 
km/s speed compared to the speed of light, we consider the Earth and Sun to be at rest in 
the CMB frame. 

To determine A, we need to specify an initial reference frame that transforms to the 
CMB frame with A. We assume that the initial frame is the particle's rest frame and that 
the transformation A takes the rest momentum { , m} to the trajectory momentum p^ = 
{"p*,p*}, i.e. A is the transformation denoted L in ([5]). We write A = L in the form 

AUp) = LUp) =SI + {1 + ^r'm'Yp' , 
A^ = 4 = Lf = m-V and A^ = L^ = 7 = m-y , (15) 

where i,k E {1,2,3} = {x, y, z} and m is the mass of the particle. 

The last thing to choose is A, the Lorentz transformation taking the quantum field's 
frame with coordinates x to the coordinates y for the particle state, x = Xy. Remarkably, A 
= 1 doesn't work. In order for the 3-momentum p to increase when the energy p* increases, 
one can choose A to be the time inversion A = diag{l, 1, 1, —1}. 

With M, A, and A determined, one obtains an expression for p. By ( TTTj) . (iT4l) . and ( TT5l) 
one finds that 

p^' = pt'(l + 40 - 47^0) ^ p^'{l - 47V) , (16) 

where is small, i.e. | |<C 1, p'^ indicates the spatial part of the effective four- momentum, 
and k G {1, 2, 3}. Also, p^ = p^ = E + mc^ is the effective total energy of the particle state, 
including the rest energy. 

To transform spectra from one gravitational potential to another, we need the time 
component of the momentum in flT6|) . One finds that 

^^ {Eo + mc^f 
E = Eo- 40 ^-^ , (17) 
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where the total particle state energy is E + mc^ = p* and Eo is the kinetic energy at </) = 0. 
The term mc^ has been cancelled from both sides of the equation. The cubic term gets one 
p*^ = Eq + mc^ from the factor p preceding the parentheses in flT4|) and flT6|) while the other 
two come from 7 = p^ /m. 

For the energy cubed term to be important, the energy Eq dX cf) = Q must be much larger 
than the rest energy mc^, since the potential is small in weak gravitational fields. In that 
case, ( TT7|) reduces to ([T]). 

It may be that some massive particles have M, A, and A determined as above to obtain 
(TT7|) . while other particles are described with different assumptions for M, A, and A. As 
mentioned in the Introduction, accelerator experiments have shown that electrons do not 
obey (1171) at least with = 0^. It may be that the comoving coordinate system for electrons 
occurs at so small a scale that the potential due to the Sun and Earth vanishes, 00 = 0, just 
as the potential due to the Galaxy and larger structures vanishes, 0oo = 0, in a comoving 
coordinate system in free fall with the Solar System. Or it could be that A is the identity and 
Mt in flTTl) is Mt = — 1 — with Mr and A as above. These choices for M, A, and A would 
give the same curved spacetime metric ( ITOl) as above, but ( fT4l) now yields p^ = p^ + 0{(l?) 
instead of ( TT7|) . meaning the energy of the particle state would not depend linearly on the 
gravitational potential 0. 

In this paper only protons are considered as cosmic ray primaries and protons are assumed 
to have potential-dependent particle state energies as obtained above for (TT71) . 

By the definition of 'weak field', as an energetic particle Eq » mc^ moves in a weak 
field, its trajectory is accelerated by only a small amount. Thus the energy Eq is very nearly 
equal to the kinetic energy all along the particle's path in weak gravitational fields. Hence, 
Eq approximates the 'trajectory energy'. 

Consider a particle entering a nonzero gravitational potential < 0. Its particle state 
energy E increases by ([1]), while its trajectory energy increases by a very small amount. 
Such a particle then has two kinds of energy: the particle state energy and the trajectory 
energy. And these can differ significantly at sufficiently high energies. Measurements of the 
particle's trajectory by time-of- flight, say, should produce the energy Eq, but, in a particle- 
state-changing collision, the particle should deposit its particle state energy E. 

3 CRs at Earth to VHE Gamma Rays at Earth via the 

SNR 

Given the observed cosmic ray spectrum and the CR's energy dependence on gravitational 
potential in ([T]), one can deduce the CR spectrum at a typical SNR. Applying astroparticle 
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physics to the colhsions of CRs and ambient matter, one can deduce the gamma ray spectrum 
produced by the colhsions at the SNR. The VHE gamma ray spectrum at the SNR and the 
assumption that gamma rays also obey a modified version of ([1]) allows one to deduce the 
expected VHE gamma ray spectrum at Earth. The process is the topic of this section. 

The observed CR fiux dNp/dEp^ (Afanasiev (Yakutsk) 1996, AGASA 2003, Auger 2007, 
Bird (HiRes) 1994, Cronin et al. 1997, Grigorov et al. (Proton) 1991, Lawrence et al. (Hav- 
erah) 1991, Nagano et al. (Akeno 1) 1984, Nagano et al. (Akeno 2) 1992, Seo et al. (LEAP) 
1991) is sectioned by energy and each section is fit. The highest energy CRs observed are 
not relevant here, so the data has not been updated for the latest values. Computer soft- 
ware (Mathematica 6.0) is used to determine best fit values of {oj, 6j, q} in log {dNp/dEp^) = 
ai + bi log (Epe) + Ci[log (£^Pe)]^ for five intervals of log (Ep^) : {log (^Pffi)i„itiai> log (^P®)finai} 
= {8.33,10.64}, {10.64,14.92}, {14.92,16.68}, {16.68,18.54}, {18.54,20.5}, where Ep^ has 
units of eV. A collection of points {log Ep^, log dNp/ dEp^}, calculated every 6 log Ep^ = 
0.1 along the five fits, is made from log{Ep^) = 8.0 to 22.0. Minor adjustments are made 
to smooth the transition from one section to the next. A third-order interpolation is con- 
structed based on the collection of points. The data and interpolating function are shown 
in Fig. m 

Eq. flT7|) from Sec. [2] relates the CR energy Ep^ at Earth (©) and the CR energy Eps 
at the supernova remnant (S). However, the expressions are messy, so, for the purpose of 
displaying the intermediate results in this section, the simpler and only slightly less accurate 
formula ([1]) is used. 

The fiux at the SNR, dNp/dEps, is related to the fiux at earth by 

(18) 

where the CR energy at null potential EpQ gives both the CR energy at Earth, Ep^ = 
Epo — 4:(I)qEpI/ TTipC'^ , and the energy of the same CR at the SNR, Eps = EpQ — AipsEpl/mpC'^. 
Many SNRs contribute to the observed CR spectrum at Earth, so the SNR is called a 'typical' 
SNR. Fig. [2] shows the CR fiux at the typical SNR using various values of (ps- The absolute 
normalization is meaningless because it depends on the unknown number of contributing 
SNRs and the unknown distances to the SNRs. 

It is assumed that VHE gamma rays are produced when CR protons strike ambient 
protons at the SNR. The VHE gamma ray fiux at the SNR, dN^/dE'js, is determined from 
the CR fiux at the SNR, dNp/dEps, using an analytical expression for secondary gamma-ray 
spectra from inelastic proton-proton interactions presented by Kelner et al. (2006). One has 

d'^N^ d^N^ dNrniiL,in„ dN„ dx 

(19) 



dNp 


7 Tvr c'-^P® 
"^Vp dEpo 


dNp 1- 


- 120e ^ 


dEps 


dEp^'^ 


dEp^ 1- 


- 1205 ^ 



V7 _ a-lM^ a^"^ collision «JVp 



dE-fsdEps dxdN collision dNp dEps dEjs 
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where d? N^ / dxdN collision = Fyix, Ep) indicates the number of gamma rays in the interval 
(x, X + dx) per colhsion and x = E'js/Eps. Also, dN collision /dNp = cnuCTinei where c is the 
speed of light, nu density of ambient protons, and ainei is the inelastic p-p cross section. It 
follows that 

dN^ /•- /E7s p \ ij, . dNp dEps 



dEjs JE-ys \Eps J dEps Eps 

where the absolute normalization A^ is selected to simplify the expression by adjusting the 
units and canceling the factor cuh, A^cuh = 1- The CR spectrum at the SNR, dNp/dEps, 
is given by (fT8|) and expressions for F^ and ainei are given by Kelner et al. (2006). The 
approach by Kelner et al. (2006) introduces a lower limit of 10^^ eV, so E'js > 0.1 TeV. The 
VHE gamma ray spectra for various values of the gravitational potential at the SNR (ps are 
displayed in Fig. HI 

A VHE gamma ray can interact with a variety of particles, low energy photons or low 
energy massive particles, to produce a massive particle such as a neutral pion or eta, an 
electron, a positron, etc. And the particle produced has practically the same energy as the 
VHE gamma ray. Even if this never happens, one can imagine it could happen anywhere 
along the gamma ray's trajectory. It is inferred that the energy of the gamma ray would 
have the same dependence on gravitational potential as any massive particle it could produce. 
Since it could produce many different massive particles, let us try using the massive particle 
Eq. ([1]) with the mass ?7i in ([1]) allowed to be some average mass, an average of the masses 
of a neutral pion or eta, electron, positron, etc. 

In effect, ([I]) is applied to VHE gamma rays with the mass m —>■ rriray allowed to vary 
freely. A best-fit value for rriray is obtained by comparing the observed VHE spectra with 
the VHE spectrum predicted for a typical SNR. None of this implies that a gamma ray has 
mass or that the imagined interactions actually occur. 

Then applying ([1]) with m —>■ ruray to the VHE gamma ray flux dN.y/dE'js at a typical 
SNR in fl2Ul) yields the expected VHE gamma ray flux dN^/dEj^ at Earth received from a 
typical SNR. One flnds 

dN, _ dN, §^ _ dN, ^-^2<Ps^ 
dEle dE^s'^ dE^sl-l2<P^-M_^ 

where the energy at null potential E'^q is a parameter that gives both the gamma ray energy 
at the Earth E'~^^ = E'jo — 4:(j)^E'j^/m'^^yC^ and the energy of the same gamma ray at the 
SNR E'js = E'jo — ^(psE'-yl/m'l^yC^. Figs. [5] and [6] gives the VHE gamma ray spectra at 
Earth for various values of (ps and rriray 

The process described in this section produces VHE gamma ray fluxes and energies from 
the observed CR spectrum at Earth given the two parameters (ps and rriray In the next 
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section the values of (f)s and rriray are found that give a best fit to the observed VHE gamma 
ray fluxes. 

4 Fitting the Observed VHE Gamma Ray Spectra 

Galactic CR protons are deflected as they travel in the Milky Way and directional information 
is lost, whereas a gamma ray trajectory aligns with the other electromagnetic emissions from 
its source. One can see where gamma rays come from. Since the CR spectrum at Earth 
represents many sources, it is appropriate to compare the VHE gamma ray spectra found in 
Sec. 3 with a VHE gamma ray spectrum obtained by combining many spectra. 

The VHE gamma ray spectra found in Sec. 3 depend on the SNR gravitational potential 

05 and the mass ruray needed to apply ([T]) to gamma rays. Given a value for (ps and and a 
value for rriray the process in Sec. 3 determines a proposed typical VHE gamma ray spectrum 
based on the CR spectrum observed on Earth. Call the VHE gamma ray spectrum obtained 
the 'proposed spectrum'; it is a function of (ps and rriray 

Since VHE gamma rays can also be produced in the lepton scenario, it is important 
to select spectra produced by the hadron scenario exclusively, with little contamination by 
VHE gamma rays produced by the lepton scenario. Among the qualifying sources are the 
spectra from the SNRs RX J1713. 7-3946 {h = 28 data points) and RX J0852.0-4622 (/i2 = 
14 data points) and from the Galactic Center ridge {hs = 9 data points) (Aharonian et al. 
2007a, 2007b, 2006). In this section, these 51 data points are combined and compared with 
the VHE gamma ray spectra found in Sec. 3. 

To compare spectra, the absolute normalization of the proposed spectrum is not changed 
while the absolute normalizations of the observed VHE gamma spectra are adjusted to 
minimize the separation of the experimental flux points to the proposed spectrum. Let 
Flux{E^,(f)s, rriray) = dN^/dE'j^Q be the proposed flux and let Flux^j{Eij) = {dN^/dE^)ij 
be the jth experimental flux measured at an energy of Eij from the ith source. Then, for 
each of the three sources, the value of Aj, z G {1, 2, 3}, is found that minimizes the average 
separation (fi(Aj) given by 



5i{Ai) = ( ;^ E [^i + log {Flux^ijiEij)) - log {Flux{Eij, (f)s, rriray))] ' , (22) 



where hi is the number of data points for the ith source. The normalization- adjusted flux 
at an energy of Eij is given by 

log (FluXijiE,,)) = Ai + log [Flux^iEij)) . (23) 
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The collection of the three adjusted spectra are now to be compared with the proposed 
spectra from Sec. 3. 

Once the absolute normalizations of the three sources are adjusted, the goodness-of-fit of 
a proposed VHE gamma ray spectrum is calculated. The goodness-of-fit function g{(f)s, ^^ray) 
is the average separation of the proposed spectrum from the Y.K — ^ = 50 observed flux 
data points. The rejected data point is the flux with a gamma ray energy of 169.79 TeV in 
the spectrum of RX J1713. 7-3946, rejected because it has low significance (Aharonian et al. 
2007a). The goodness-of-fit function g{4'Si'^Tay) is defined by 

I 1 2\ "' 

9{4>s, rriray) = ^ _ ^ Y. [log {FluXij{Eij)) - log {Flux{Eij, (ps, mray))f ■ (24) 

By varying (ps and rriray and calculating g^cps^i^ray) for each, a best fit is found. The best 
fit and some close runnersup are given in Table [H 

The best fit has a goodness-of-fit of 0.074 with (ps = 0.110®, rriray = 0.45m^o, where 
m^o is the mass of a neutral pion and 0® = —1.06 x 10~^ is the Newtonian potential at 
Earth due to the Sun and Earth referenced to a null potential at great distances from the 
Solar System, see (TT2l) . The best fit proposed VHE gamma ray spectrum and the three 
normalization- adjusted observed spectra including the point at 169.79 TeV are plotted in 
Fig. El 

Tables [2], [3] and |4] contain sample numerical values for the various quantities discussed in 
this article. Wherever needed, the tables use the best fit value of the gravitational potential 
at the SNR, 05, and the best fit value for rriray 

5 Discussion and Conclusion 

The explanation of the VHE gamma ray spectra of typical SNRs presented in this paper is 
based on the dependence of particle state energy on gravitational potential as displayed in 
([1]). A proton or a gamma ray in the near-null potential at an SNR has less energy than it 
has in the gravitational potential at the Earth. The quality of the fit to the data in Fig. [7] 
supports this explanation. 

The quality of the fit supports the application of formula ([T]) to photons even though 
photons are massless and ([T]) applies to massive particles. It is argued that the photon 
should have particle states with energies almost the same as the energies of various massive 
particle states. The coincidence of photon particle state energy with massive particle state 
energy is supported by the VHE gamma ray spectra. 

For a given source of VHE gamma rays, the distance to the source and the number of 
CRs each contributes to the CR spectrum at Earth are not knowable to sufficient accuracy 
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to be useful in setting the absolute normalization of the VHE gamma ray spectra. One might 
include as free parameters the AjS needed to adjust the absolute normalizations of the three 
VHE gamma ray spectra. 

One of the two parameters needed to obtain an expected VHE gamma ray spectrum from 
the observed CR spectrum at Earth is the mass rriray Since the best fit mass rriray is less 
than the neutral pion mass, there must be some interaction that allows a VHE gamma ray 
to transfer almost all its energy to a particle with less mass than the neutral pion. Keep in 
mind that transferring energy to electrons and positrons does not contribute to rriray because 
they do not obey ([1]). 

The best fit gravitational potential at the SNR is about ten percent of the potential at 
the Earth's surface. Using this to imply something about the density of, the dimensions 
of, and the distance to the clouds of ambient matter that produce the VHE gamma rays is 
beyond the scope of this article. 

Consider Fig. 7. Draw a straight line on the graph from a point near 0.3 TeV to a point 
near 30 TeV, on the graph this is from logEj^) = 11.5 to 13.5. Note the 'bump' in the VHE 
gamma ray spectra. The same bump is seen in all three observed spectra which implies that 
some common mechanism is in place to produce the bump. That the same bump can be 
replicated by the best fit values of (ps and rriray shows how well the explanation relates the 
CR spectrum at Earth with the VHE gamma ray spectrum at Earth. 

The tightness of the fit in Fig. [7] suggests that the feature found in the VHE gamma ray 
spectrum of SNR RX J1713. 7-3946 and the other SNRs plotted in Fig. [7] exhibits particle 
behavior consistent with the simple expression ([T]) of a gravitational potential-dependent 
relationship between CR and VHE gamma ray energies at Earth and the energies at the 
SNR. 
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Table 1: There are low goodness-of-fit values g{4'S)iT^ray) for proposed VHE gamma ray 
spectra with rriray from 0.35 to 0.55m^o and with ^^ from 0.09 to 0.1300, where rriray is the 
average mass of possible gamma ray produced particles and 05 is the gravitational potential 
at the SNR. The constant 00 is the potential at the Earth's surface due to the Earth and 
Sun referenced to zero potential at infinite distance, see (IT2|) . The As adjust the absolute 
normalizations of the three observed spectra to each proposed spectrum. 



# 


^(0s,m™y) 


rriray/m^o 


0s/0e 


Ai 


A2 


A3 


1 


0.074 


0.45 


0.11 


1.75 


1.62 


0.91 


2 


0.076 


0.35 


0.11 


1.67 


1.54 


0.83 


3 


0.077 


0.45 


0.09 


1.72 


1.60 


0.89 


4 


0.078 


0.55 


0.11 


1.81 


1.69 


0.97 


5 


0.078 


0.45 


0.13 


1.78 


1.65 


0.93 


6 


0.079 


0.55 


0.09 


1.79 


1.67 


0.96 


7 


0.080 


0.35 


0.09 


1.63 


1.50 


0.80 


8 


0.080 


0.35 


0.13 


1.71 


1.57 


0.85 


9 


0.082 


0.55 


0.13 


1.83 


1.71 


0.99 
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Table 2: The particle state energies of 11 CR protons in the gravitational potential of the 
Earth (©), in a null potential (0), and in the potential at an SNR (5), providing examples 
of Eq. ([1]) which show how particle state energy depends on gravitational potential. For the 
potential at the SNR, (ps, the best fit value from Tabled] is used. 



# 


Epe (TeV) 


Epo (TeV) 


Eps (TeV) 


1. 


1.00 


0.96 


0.96 


2. 


2.00 


1.74 


1.77 


3. 


4.00 


2.87 


2.99 


4. 


8.00 


4.27 


4.68 


5. 


16.0 


5.94 


7.05 


6. 


32.0 


7.94 


10.6 


7. 


64.0 


10.4 


16.3 


8. 


128. 


13.4 


26.0 


9. 


256. 


17.1 


43.4 


10. 


512. 


21.7 


75.6 


11. 


1024. 


27.5 


137. 
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Table 3: The energies of 11 VHE gamma rays in the gravitational potential of the Earth 
(©), in a null potential (0), and in the potential at an SNR (S). While ([T]) only applies to 
massive particles, it is applied here for massless gamma rays by considering virtual though 
possible processes that produce massive particles with essentially the same energies as the 
gamma rays. The energy dependence of a gamma ray should average the energy dependences 
of these virtual particles. There are many possible processes, i.e. inverse pion decay, the 



Compton effect, pair production, etc. An average mass called rriray is used in ([T]) and is 
varied to fit the observed data. The best fit values of rriray and (ps in Table [1] are used with 
([1]) to obtain the values of E'Jq and -^7© in this table. 



# 


E^s (TeV) 


Ejo (TeV) 


E7e (TeV) 


1. 


0.164 


0.159 


0.204 


2. 


0.301 


0.275 


0.514 


3. 


0.509 


0.417 


1.25 


4. 


0.795 


0.566 


2.65 


5. 


1.20 


0.722 


5.04 


6. 


1.80 


0.895 


9.12 


7. 


2.77 


1.10 


16.3 


8. 


4.42 


1.35 


29.3 


9. 


7.35 


1.66 


53.6 


10. 


12.9 


2.05 


100. 


11. 


23.3 


2.54 


191. 
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Table 4: The CR and VHE gamma ray fluxes at Earth and at the typical SNR for the 
CRs and gamma rays in Tables [2] and [31 The CR flux at Earth, dNp/dEp^, represents 
experimental data. The CR flux at the SNR follows from the CR flux at Earth and particle 
state energy dependence, flTSl) . The gamma ray flux at the SNR is the result of integrating 
over the CR proton energy at the SNR using the method of Kelner et al. (2006). The gamma 
ray flux at Earth is found by (^11) . The best fit values of rriray and (f)s in Table [T] are used to 
obtain the calculated results in columns 3, 4 and 5. 



# 


dNp/dEp^ 


dNp/dEps 


dN-f/dE-fs 


dN-f/dE-f^ 


1 


2.18 X 10-13 


2.43 X 10-13 


2.87 X 10-11 


1.69 X 10-11 


2 


3.48 X 10-14 


4.78 X 10-14 


5.23 X 10-12 


1.88 X 10-12 


3 


5.56 X 10-15 


1.08 X 10-14 


1.15 X 10-12 


2.73 X 10-13 


4 


8.87 X 10-1^ 


2.49 X 10-15 


2.92 X 10-13 


5.37 X 10-14 


5 


1.41 X 10-16 


5.52 X 10-16 


7.70 X 10-14 


1.21 X 10-14 


6 


2.25 X 10-1' 


1.14 X 10-16 


1.90 X 10-14 


2.69 X 10-15 


7 


3.58 X 10-is 


2.19 X 10-1^ 


4.11 X 10-15 


5.38 X 10-16 


8 


5.70 X 10-19 


3.98 X 10-1^ 


7.70 X 10-16 


9.55 X 10-1^ 


9 


9.06 X 10-20 


6.93 X 10-19 


1.33 X 10-16 


1.58 X 10-1^ 


10 


1.44 X 10-20 


1.17 X 10-19 


2.23 X 10-1^ 


2.59 X 10-1^ 


11 


2.27 X 10-21 


1.92 X 10-20 


3.70 X 10-is 


4.22 X 10-19 
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A Problems 

[Answers can be found in this section in the LaTeX source document hidden by %s.] 

1. (a) Use the values of Epo hsted in Table [2] to calculate £'p© and Eps using formula ([T]). 
Find the percent differences with the values of Ep^^ and Eps found using (TTTj) . (b) Similarly, 
for the i?7o listed in Table [3l find the percent differences of the values of Ej^ and Ejs found 
by (P and (^. 

2. Wherever needed, use the best fit values from Tabled] in this problem, (a) Use the energies 
listed in Table [21 the fluxes in Column 2 of Table H] and (TT8|) to verify the CR fluxes at the 
SNR, dNp/dEps in column 3 of Table HI (b) Likewise, use (12T|) to verify the dN'j/dE'yt^ 
listed in Column 5 of Table HI 

3. The 11 gamma ray energies E'ys in Table [3] are each 17% of the 11 CR energies Eps 

in Table [2] in keeping with the ^-functional approach that assumes secondary gamma ray 

energies are a fixed fraction, k, = 0.17, of incident CR proton energies. Given a proton flux 

at the SNR, dNp/dEps, from Table H] and the formula for the inelastic cross section atnei, 

Eq. (79) from Kelner et al.[Tl], with the normalization chosen to cancel the factors c and 

hh, one can obtain gamma ray fluxes dNj/dE'js using the (5- functional approach. One finds 

that 

dNp 



dN^ 
dE^ 



AcnHH(Tinel{Eps, ,^ 

dEps 



'J &— functional 

where the normalization can be adjusted to match the values in Table H] (I use Acnu = 20.) 
Compare the resulting values with the ones listed in Table H] which are found by integration, 
see Eq. (71) from Kelner et al.p7 



4. The observed CR spectrum in the region of interest between the dashed lines in Fig. 1 
is approximately a power law, dNp/dEp^ = Ap^Ep"^^^, where Ap,^ is a constant, (a) Find 
the range of the parameter Epo in the region of interest, (b) Using Epo as a parameter, 
or otherwise, find an expression for the CR spectrum at the SNR, dNp/dEps- (c) Use the 
5— functional method, described briefly in Problem 3, to obtain an expression for the VHE 
gamma ray spectrum at the SNR, dN'y/dE'js- [It may be useful to use E'jo as a parameter 
and find the range of E'Jq that corresponds to the range of E'js (= i^Eps).] (c) Find an 
expression for the VHE gamma ray spectrum at the Earth. How does this result compare 
with the result using the analysis of Kelner (2006) displayed in Fig. 7? 

5. Can you give an explanation of why the best fit mass niray is about half of a neutral pion 
mass? 



A PROBLEMS 



24 



10 11 12 13 14 15 16 17 18 19 20 21 



> 



3 

O 



-5 



-10 



-15 



-20 



-a 

Z; 

^ -25 



-30 



-35 



-40 



Interpolation 
(Ojfset by -2) 



) Knee 




• AGASA 

• Akeno 1 
n Akeno 2 
V Auger 
■ Haverah 



T HiRes 
A LEAP 
O Proton 
o Yakutsk 



- -5 



- -10 



- -15 



- -20 



- -25 



- -30 



- -35 



- -40 



8 9 10 11 12 13 14 15 16 17 18 19 20 21 
log(Ep ) with Ep in eV 



Figure 1: The Cosmic Ray Spectrum at Earth. The dashed hues indicate the region of 
the observed cosmic ray (CR) spectrum that corresponds to the observed VHE gamma ray 
spectrum. In that region the CR spectrum is a near-pure power law without any distinctive 
feature. The 'Knee' in the CR spectrum, where the spectral index changes, occurs at a 
higher energy and should have little effect on the observed VHE gamma ray spectrum. 
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Figure 2: The Expected Cosmic Ray Spectrum at a Typical SNR. The spectrum bends in the 
region corresponding to the VHE gamma ray spectrum. Compare this with the same region 
of the CR spectrum at Earth, see Fig. [H The bend in the CR flux at the SNR reflects the 
CR proton energy dependence on gravitational potential, ( TTSj) . The displayed curves are for 
potentials (ps differing by successive factors of 2. By formula ([1]), the CR energy at Earth is 
Ep^ = EpQ — 4:(j)^Epl/mp'^ and at the SNR, Eps = EpQ — 4:(j)sEpl/mp'^, so at high energies, 
we have Eps = Ep^(f)s/4><s- It follows that, at high energies, a given CR flux appears at Eps 
values separated by about A{\ogEps) = log 2 = 0.30. 
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Figure 3: Contributions to the VHE Gamma Ray Spectrum from the p-p Cross Section 
and the Energy Transitions of the CR and Gamma Ray Fluxes . [A modified 5-functional 
approach is used for this rough sketch: instead of E'js = 0.17i?p5, we employ E'Jq = 0.17EpQ.] 
(a) log ainei for rriray = O.Sm^o. The log of the inelastic p-p cross section slopes steadily upward 
in the region corresponding to the observed VHE gamma ray spectra, (b) log (dEp^/dEps) 
+ log (dE'js/dE'y^) for (ps = 0.10®. Relating the CR flux at Earth to the flux at the SNR 
and relating the VHE gamma ray flux at the SNR to the flux at Earth introduces the factors 
(1 - 120e 5|)/(1 - 1205 51) and (1 - 120^ J%)/(1 - 120e ^), respectively. See 

p p ray ray 

flTSjl and fl2Tl) . At both high and low energies the product of the factors is unity, but in the 
region of the VHE gamma ray spectrum the log of the product of the factors changes by up 
to 0.5 (= log 3) depending on rriray and 05. 
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Figure 4: The Expected VHE Gamma Ray Spectrum, at a Typical SNR. CR protons must 
have more energy than the secondary gamma rays they produce. So the gamma ray spectra 
are calculated by integrating over those CR protons that have more energy. The calculation 
follows the method of Kelner et al. (2006). Note how the gamma ray spectrum at the SNR 
in this graph mirrors the CR spectrum at the SNR in Fig. [2], except that the gamma ray 
energies are lower than the CR energies. 
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Figure 5: The Expected VHE Gamma Ray Spectrum, at Earth; Dependence on the SNR 
Potential (ps- At high energies, just as expected from Fig. [3)3, taking the CR spectrum at 
Earth backwards to the SNR cancels taking the VHE gamma ray spectrum at the SNR 
forward to Earth, as far as dependence on (ps is concerned. The only dependence on (ps is 
in the region of interest, the region of the feature in the observed VHE gamma ray spectra. 
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Figure 6: The Expected VHE Gamma Ray Spectrum at Earth; Dependence on Mass rriray 
Just as with the dependence on the SNR potential (j)s, the dependence on the quantity ruray 
is confined to the region of interest. The reason is the same, in the hmit of large or small 
values of rriray, the factor whose log is plotted in Fig. [3)d approaches unity. The curves are 
drawn left to right in the order that the m^ay^ are listed. 
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Figure 7: T/ie i?esi Fit to the Observed VHE Gamma Ray Spectrum, at Earth. The graph 
shows the proposed VHE gamma ray spectrum at Earth evaluated in the region of interest 
for the best fit values in Tabled! The data points are compiled from the published spectra of 
SNRs RX J1713.7-3946 and RX J0852.0-4622 and from the Galactic Center ridge (Aharonian 
et al. 2007a, 2007b, 2006). The absolute normalizations of the three spectra are adjusted 
by the values Aj listed in Table [U see fl2^ . The three highest energy data points are 
outnumbered by the 48 data points of lower energy, so fitting the three highest energy data 
points so closely is remarkable. 



